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Artificial light is globally one of the most widely distributed forms of anthropogenic pollution. However, 
while both the nature and ecological effects of direct artificial lighting are increasingly well documented, 
those of artificial sky glow have received little attention. We investigated how city lights alter natural regimes 
of lunar sky brightness using a novel ten month time series of measurements recorded across a gradient of 
increasing light pollution. In the city, artificial lights increased sky brightness to levels six times above those 
recorded in rural locations, nine and twenty kilometers away. Artificial lighting masked natural monthly 
and seasonal regimes of lunar sky brightness in the city, and increased the number and annual regime of full 
moon equivalent hours available to organisms during the night. The changes have potentially profound 
ecological consequences. 



The widespread use of artificial lighting to illuminate the night-time environment has raised significant 
concerns regarding its aesthetic^"^, environmentaP"^, and human health impacts^'^. Yet despite measures 
to improve economic efficiency and mounting pressure from the astronomical community, light pollution 
persists, continues markedly to grow in spatial extent^, and is broadening in its spectral characteristics^. 

Aesthetically, sky glow (artificial light which is scattered and reflected back to earth by the atmosphere) 
obscures our view of natural starlight and moonlight^°'^\ In 2001 it was estimated that 60% of the world's 
population inhabited areas affected by artificial sky glow^ leading astronomers to campaign for darker skies in 
urban and suburban regions ^ and the establishment of the first international dark sky parks for recreational 
stargazing^^. Associated with increasingly urbanised lifestyles, this loss of naturally lit sky arguably contributes to 
the 'extinction of experience' or the disengagement of people from the natural environment^^, compounding 
efforts to engage the public on issues of environmental concern. 

In addition, sky glow may be having environmental impacts that have so far remained largely unexplored, while 
those of direct artificial lighting are much better understood and increasingly well documented. The impacts of 
direct artificial light include navigational disorientation of birds^^'^^, turtles^^'^^ and invertebrates^^'^^, the displace- 
ment and concentration of selected foraging habitat in bats^°'^^ and birds^^, the extension of foraging times in 
birds^^"^^, and changes to the composition of invertebrate communities^^. Added to this, the potential effects of sky 
glow include artificially increasing the amount of ambient light and time available for foraging^^'^^, interfering 
with celestial light patterns used in navigation^^, and the disruption of temporal activity patterns linked to daily, 
monthly and annual lunar rhythms in sky brightness^^"^^. Understanding how sky glow alters natural regimes of 
sky brightness is therefore fundamental to investigating its environmental impacts. 

Natural regimes of night-time sky brightness are influenced by the harmonic movements of the astronomical 
bodies. During each month, sky brightness resulting from lunar rhythms varies from its highest when the moon is 
full and above the horizon, to its lowest when the moon is new^^ At seasonal time scales, variation in the altitude of 
the full moon results in winter nights being more brightly lit by the lunar disc than summer nights, while at longer 
time scales the 18.6 year cycle of variation in the moon's orbit also affects night-time sky brightness regimes. 
These natural predictable regimes in sky brightness provide organisms with a lunar clock with which to syn- 
chronise activities. In addition, the variation in night-time length at higher latitudes gives rise to predictable 
annual regimes in the amount of moonlit hours available to conduct activities such as foraging and reproduction 
in any given month. Artificial light pollution may interfere with monthly regimes in sky brightness, and annual 
regimes in the amount of sufficiently lit foraging time available to nocturnal species. 

Here we determine how artificial light pollution alters natural sky brightness regimes. First, we assess whether 
artificial lighting changes monthly and seasonal regimes in sky brightness compared to those which occur 
naturally as a result of lunar rhythms. Second, we determine how this decoupling affects the number of full 
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moon brightness equivalent hours available for nocturnal activity 
across the year. Throughout, we are considering only the diffuse light 
scattered across the sky dome, rather than the direct light from the 
lunar disc itself 

As a case study, sky brightness was monitored for 10 months at 
three locations which varied in proximity to Plymouth city centre, 
UK (Figure 1). Changes in sky brightness resulting from temporal 
variation in moon altitude and phase were then modelled for the 
same period. The observed values and model predictions of sky 
brightness were compared to understand how artificial light pol- 
lution alters monthly and annual regimes of night-time sky bright- 
ness at each of the three locations. The number of full moon 
equivalent sky brightness hours available for activity per night in 
each month was compared between the three sampling locations. 
This was quantified as the number of hours when the sky was 
brighter than the 3""^ quartile observation made when the moon is 
between three-quarters and full (<45° phase angle) at a location 
where sky brightness was predominantly governed by celestial 
bodies. 

Results 

Night time sky brightness during the observation period was on 
average 8 times greater in Plymouth city centre ( 0.048 ± 0.001 SE 
cd m-^ SD = 0.051) compared to locations 8.6 (0.006 ± 0.003 SE cd 
m-^ SD = 0.012) and 19.7 (0.006 ± 0.0005 SE cd m'^ SD = 0.022) 
km away (SE = Standard Error, SD = Standard Deviation). Monthly 
variation in sky brightness due to the orbit of the moon around the 
earth (phase) and seasonal variation in maximum sky brightness due 
to variation in the altitude of full moon (Figure 2d) were easily 
recognised in both the 19.7 km (Figure 2a) and 8.6 km (Figure 2b) 
sampling locations. By contrast, natural rhythms of lunar sky bright- 
ness were not evident in observations recorded from the sky over 
Plymouth city centre (Figure 2c). Observations of sky brightness on 



nights of lunar activity (when the moon was above the horizon) over 
the sampling period were significantly correlated with predicted 
lunar sky brightness at all three locations (Pearson correlation; city 
centre:r = 0.16,df = 1186, p < 0.0001; 8.6 km: r = 0.70, df = 1106, 
p < 0.0001; 19.7 km: r = 0.59, df = 1036, p < 0.0001) (df = degrees 
of freedom). However, the strength of this correlation was markedly 
lower in Plymouth city centre compared to both the 8.6 km and the 
19.7 km sampling locations, indicating that observed sky brightness 
regimes were less synchronised with predicted regimes in lunar sky 
brightness in the city (Figures 2 and 3). Indeed, while the two rural 
locations displayed a clear linear increase of observed with predicted 
sky brightness (Figure 3a and b, note logarithimic y axis), observed 
scaled with predicted sky brightness in a polygonal relationship in 
Plymouth city centre (Figure 3d). Here, minimum observed sky 
brightness was determined by predicted lunar sky brightness 
regimes, while maximum brightness was not. The variation in 
observed sky brightness above the minimum is likely a result of the 
influence of artificial lights and their interaction with cloud which 
can amplify sky brightness in the city by reflecting artificial light back 
towards the earth's surface^^. 

The location of this study in northern temperate latitudes means 
that astronomical night length varies substantially throughout the 
year (Figure 4). Hence the average nightly number of full moon 
equivalent hours was higher in winter across all three locations. 
The difference in the sky brightness regimes in Plymouth city centre 
resulted in significantly more full moon equivalent hours per night, 
compared to the two rural locations (ZAP, hours —site: df = 4,%^ = 
535.7, p < 0.0001) with the magnitude of this increase being more 
pronounced in winter than summer months (Figure 4). In the city, 
the presence of artificial lighting means that the annual regime in the 
average nightly hours of full moon equivalent sky brightness more 
closely approximates the regime in astronomical night length, while 
in the countryside this remains linked to regimes of lunar periodicity 




Figure 1 | Sampling locations in relation to artificial light pollution. SQM-LU-DLs were deployed in Plymouth city centre (A), 8.6 km from the city 
centre (B), and 19.7 km from the city centre (C). Image acquired by the Suomi National Polar-orbiting Partner satellite on March 27, 2012. Downloaded 
from http://earthobservatory.nasa.gov/IOTD/view.php?id= 78674 on 13* December 2012, georectified using ArcGIS v9.2. 
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Figure 2 | City lights disrupt regimes in lunar night- time sky brightness. Observed regimes of night-time sky brightness (grey points) recorded 19.7 km 
(a), 8.6 km (b) and 0 km (c) from Plymouth city centre are compared to the predicted natural regime of lunar sky brightness (d). Points are plotted with 
transparency to provide information about density (darker regions correspond to denser aggregations of points). Observations are plotted for the period 
between 01/10/201 1 and 31/07/2012. Note that the gap in the data around 2456100 julian date occurs because, at our locations, astronomical night does 
not occur during June. Each point represents an hourly recording of sky brightness (a-c) or predicted relative sky brightness (d), with recordings made 
when the moon was below the horizon, and when SQMs threshold at maximum brightness omitted (see methods). 



(Figure 4). This difference in the regime of average nightly hours of 
full moon equivalent sky brightness was significant between sites at 
the 95% confidence level (ZAP, hours~site*month: df = 32, %^ = 
46.6, p = 0.05). 

Discussion 

The results of this study demonstrate that artificial sky glow from city 
lights substantially alters natural regimes of night-time sky bright- 
ness. The presence of artificial lights increases the average intensity of 
night-time sky brightness, masks natural regimes resulting from 
lunar periodicity, and increases the number of full moon equivalent 
hours of sky brightness throughout the night. This latter effect is 
greater in winter than in summer at higher latitudes, resulting in a 
change in the annual regime of full moon equivalent hours of sky 
brightness. 

The astronomical community have typically led the campaign for 
the reduction of artificial sky glow on aesthetic and scientific 
grounds\ This study has highlighted how artificial lights alter natural 



regimes of sky brightness at different temporal scales. These changes 
may be having deleterious impacts on a variety of ecological pro- 
cesses which are informed by the periodicity in the brightness of the 
night sky. Increasing the hours of full moon equivalent sky bright- 
ness per night may extend or reduce the time available to organisms 
to perform activities. For example, a very large number of species are 
nocturnal or facultatively nocturnal, with foraging activity either 
restricted to or increased on moonlit nights. The artificial light from 
nearby urban centers may extend the number of hours such species 
spend foraging at night, and this has indeed been demonstrated for 
shorebirds feeding on mudflats^^'^^. Conversely, prey species may 
experience increased predation rates or, where mobile, reduced for- 
aging times as a consequence of predator avoidance. For example, the 
diel vertical migration of zooplankton is a light intensity dependent 
behaviour which is thought to reduce predation risk^^. However this 
migration has been observed to be suppressed in freshwater lakes 
above which the night-time sky is artificially brightened^°. Many 
behaviours are often synchronized to short term variations in the 
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Figure 3 | City sky brightness correlates weakly with predictions based on 
lunar activity. The correlation between observed sky brightness and the 
predicted natural regime of lunar sky brightness is presented for 19.7 km 
(a), 8.6 km (b) and 0 km (c) from Plymouth city centre. Points are plotted 
with transparency to provide an information about density (darker regions 
correspond to denser aggregations of points). Observations are plotted for 
the period between 01/10/2011 and 31/07/2012. Each point represents an 
hourly recording of sky brightness, with recordings made when the moon 
was below the horizon, and when SQMs threshold at maximum brightness 
omitted (see methods). 

illumination of the moon's disk throughout each 30 day phase per- 
iod, a phenomenon which is widespread throughout the animal 
kingdom^^'^^'^^"^^. Hence, the short term destabilization of monthly 
lunar sky brightness regimes by artificial lights may disrupt patterns 
in foraging and mating linked to these phase periods. The increase in 
the magnitude of sky glow by artificial lights in cities is also likely to 
affect these behaviours and, in addition, may interfere with the nav- 
igation of many species which use celestial patterns of light to ori- 
entate themselves in the night-time environment^^'^^. For example, 
artificial light has been shown to disrupt the orientation and there- 
fore migration direction of birds^^, and concerns have been raised 
about the interference of artificial light pollution with celestial polar- 
ized sky light patterns^^ which many insects^^"^° and birds^^'^^ use to 
navigate. 
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Figure 4 | City lights increase the nightly number of full moon equivalent 
hours more during winter months than summer months. The monthly 
variation in the number of hours per night when sky brightness was 
observed to be greater than the 3'^'^ quartile sky brightness observed in each 
month at a dark sky location (19.7 km from the city centre) when the 
moon is close to full (<45° phase angle). Monthly average night length in 
hours is presented to illustate the natural annual variation in night length 
at the study locations. Error bars represent standard errors. Note that at the 
latitude of our locations, astronomical night does not occur during June. 

Using a 10 month time series of sky brightness measurements, this 
study has revealed how artificial light pollution can alter the natural 
regimes of lunar sky brightness which regulate a variety of beha- 
viours in many species. Yet few investigations have been conducted 
into what could be a plethora of potential ecological impacts of 
artificial sky glow. Understanding the environmental consequences 
of this already widespread^ and growing^ form of anthropogenic 
pollution is clearly essential. 

Methods 

Measuring night-time sky brightness. To achieve a gradient of light pollution, sky 
brightness was measured in Plymouth city centre, UK (50 22' 28.48" N: 
408' 19.04"W), 19.7 km from the city centre in a 'dark' rural location (50 21' 57.6"N: 
4 24' 53.01" W), and at an intermediate rural location 8.6 km from the city centre (50 
19' 01.72"N:4 13' 11.32" W) (Figure 1). Plymouth is a city covering 79.78 km^witha 
population of 240 720 people and approximately 38 000 municipal street lamps. Our 
rural sampling locations were situated on the coastline to the west of the city in order 
to minimize encroachment of artificial light from other urban centres. Sky brightness 
in magnitudes per square arc second (MPSAS) was measured every hour at each of the 
three locations using Sky Quality Meter Data Loggers (SQM-LU-DL, Unihedron) 
from the 01/10/2011 to the 31/07/2012. These values were then converted to cd m~^ 
using [cd m"^] = io.8''10^''10^"° '^''^'^^^'^^^^ which inverted the scale of the measurement 
so that it more intuitively represented sky brightness (brighter skies attain higher 
values). Each SQM -LU-DL was customised to record sky brightness directly up into 
the sky (elevation angle above the horizon of 90°) with a Half Width Half Maximum 
(HWHM) angular sensitivity of ~ 10°. Each SQM -LU-DL was switched off in order 
to download data and recharge batteries for no more than seven days on two 
occasions during the sampling period. 

Predicting night-time sky brightness. The natural regime of lunar sky brightness at 
zenith (in nL) was predicted for the time of each data point recorded using equations 
20, 21, 3, and 15 of the moon light sky brightness model developed by Krisciunas & 
Schaefe^^. The model estimates the lunar brightness of a given sky position above a 
specified geographical location at a given time under clear conditions based on the 
brightness of the moon (/) estimated as a function of its' phase angle (a), the angular 
distance between the moon and the sky position being monitored ( p), the scattering of 
light across that distance (/), and the extinction of light as it passes through the 
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atmosphere (k). The input parameters required to make the calculations are a, k, and 
the zenith distance of the moon (Zm) and the zenith distance of the sky position (Z) 
which are used to calulate p. Z is the angle between the zenith and the sky position for 
which calculations of brightness are derived. In this case the sky position observed by 
the SQM-LU-DL's was the zenith and therefore Z was always 0°. Zm was calculated as 
the moons altitude subtracted from the zenith (90°), hence p was allways equal to Zm. 
The moons' altitude was calculated from the right ascension and declination of the 
moon relative to the geocentric coordinates of each of our SQM-LU-DL locations 
using equations 8.3 and 8.4 of Meeus^^ respectively, a was estimated for each time 
point sampled using equations 31.2 and 31.3 of Meeus^^. The model accounts for the 
lunar opposition effect by multiplying sky brightness values by 1.3 when the moon is 
less than 7° phase angle^^. k is variable in both time and space, and depends on 
atmospheric conditions and height above sea level. Hence we assumed /c to be a 
constant value of 0.33, equivalent to sea level k under idealised conditions. In 
addition, the model does not account for variation in meterological conditions such as 
cloud, which can have a strong impact on sky brightness. This meant that the model 
predictions of lunar sky brightness were unlikely to match our observations in terms 
of magnitude. However, because we were interested in comparing only the correlation 
between temporal trends in observed sky brightness versus lunar sky brightness, this 
difference does not affect our conclusions. To avoid confusion, predicted lunar sky 
brightness is presented as an index between 0 and 1, with 1 corresponding to 
maximum predicted sky brightness during the observation period. 

Statistical approach. All data recorded outside of astronomical night (when the sun is 
above 18° below the horizon) were omitted to remove the influence of scattered 
sunlight on the analysis. This resulted in an absence of readings during June due to the 
lack of astronomical night at the latitude sampled when day lengths were at their 
maximum. 

The degree of sychronicity between observed and predicted natural lunar regimes 
of sky brightness was tested independently for each sampling location using Pearson 
correlations. A known problem occurs with early firmware versions on SQM-LU-DLs 
leading to recorded sky brightness thresholding, resulting in frequent recordings of 
minimum brightness which introduces a strong degree of leverage into the dataset. All 
records of maximum readings for each SQM-LU-DL were therefore omitted from 
Figures 2 and 3, and correlative analysis, to avoid this effect. The output from the 
lunar sky brightness model does not calculate sky brightness due to starlight, and 
assuming a constant value of 21.6 MPSAS on the numerous recordings during which 
the moon was below the horizon introduced a strong degree of leverage into the data. 
Therefore all recordings made when the moon was below the horizon were omitted 
from these analysis also. 

To demonstrate how changes in sky brightness regimes affect the number of hours 
of full moon equivalent sky brightness per night throughout the year, we re- analysed 
the SQM recordings including data points which had previously been omitted due to 
frequent recordings of minimum sky brightness and periods when the moon was 
below the horizon. For each sampling location, we then calculated the number of 
hours per night when sky brightness was above that when the moon is close to full and 
compared this between sampling locations. For the purposes of this analysis, the 
threshold for close to full moon sky brightness was defined as the 3''' quartile 
observation made when the moon is between three-quarters and full (<45° phase 
angle) in each individual month at the 'dark' (19.7 km) location. This was defined 
independently for each month because the brightness of the full moon varies with its 
maximum altitude throughout the year (Figure 2d). Observations made during the 
third quarter to full moon were used to define full moon equivalent brightness 
because the number of observations made when the moon is full varies between 
months, with few observations occuring in summer months. The 3""^ quartile obser- 
vation was used as the cut off in each month to ensure full moon equivalent sky 
brightness was defined by the brightest nights, while avoiding the influence of outliers 
which can have a large effect on the mean. 

The nightly number of hours where sky brightness was greater than average full 
moon sky brightness was compared between months (factor Month) and sampling 
locations (factor Site), and tested for a difference in the annual regime between 
sampling locations (Month*Site interaction). To account for the many nights that 
contained 0 hours when the sky was brighter than sky brightness during a full moon, 
we performed our analysis using Zero Adjusted Poisson Regression (ZAP). 
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